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Abstract 

We calculate the full electroweak one-loop corrections to the decay of the CP- 
odd Higgs boson A into scalar fermions in the minimal supersymmetric extension 
of the Standard Model. For this purpose many parameters of the MSSM have to 
be properly renormalized in the on-shell renormalization scheme. We have also 
included the SUSY-QCD corrections. For the decay into bottom squarks and tau 
sleptons, especially for large tan j3, the corrections can be very large making the per- 
turbation expansion unreliable. We solve this problem by an appropriate definition 
of the tree-level coupling in terms of running fermion masses and running trilinear 
couplings Af. We also discuss the decay of heavy scalar fermions into light scalar 
fermions and ^4°. We find that the corrections can be sizeable and therefore cannot 
be neglected. 



1 Introduction 



The search for a Higgs boson is the primary goal of all present and future high energy 
experiments at the TEVATRON, LHC or an e + e~ Linear Collider. Whereas the Standard 
Model (SM) predicts just one Higgs boson, with the present lower bound of its mass 
mil > 114.4 GeV (at 95% confidence level) [1], extensions of the SM allow for more Higgs 
bosons. In particular, the Minimal Supersymmetric Standard Model (MSSM) contains 
five physical Higgs bosons: two neutral CP-even (h° and H°), one neutral CP-odd (^4°), 
and two charged ones (H ) [2, 3]. The existence of a charged Higgs boson or a CP- 
odd neutral one would give clear evidence for physics beyond the SM. For a discovery 
precise predictions for its decay modes and their branching ratios are necessary. In case 
supersymmetric (SUSY) particles are not too heavy, the Higgs bosons can also decay 
into SUSY particles (neutralinos charginos xti sfermions f m ), H°,A° — > x?x° 

= 1...4), H°,A° -> xtxT (Kl = 1, 2), H°,A° - f m f n (m,n = 1, 2), H± - xUl 
H ± — > fmfn- At tree-level, these decays were studied in [4, 5]. In particular, the 
branching ratios for the decays into sfermions, H°, A — > f m f n , can be sizeable depending 
on the parameter space [6, 7]. The SUSY-QCD corrections to the decays into sfermions 
have also been calculated [8, 9]. The corrections to the decays into neutralinos [10] or 
charginos [11] due to fermion/sfermion exchanges have also been found non-negligible. 

In this paper, we study in detail the decay of the CP-odd Higgs boson A into two 
sfermions, A — > fif%. In particular, the third generation sfermions tj, 6j, and fj are 
interesting because one expects them to be lighter than the other sfermions due to their 
large Yukawa couplings and left-right mixings. Since A only couples to /l— /r (left-right 
states of /), and due to the CP nature of A , A — > fa fa vanishes. (This is valid also 
beyond tree-level for real parameters in the MSSM.) We will calculate the full electroweak 
corrections in the on-shell scheme. Owing to the fact that almost all parameters of the 
MSSM have to be renormalized in this process and hence a large number of graphs has to 
be computed, the calculation is very complex. Despite this complexity, we have performed 
the calculation in an analytic way. We have also studied the crossed channel J2 — > fiA°. 
Some important numerical results, especially for the decays of A into the third generation 
squarks and the corresponding crossed channels have already been shown in [12]. 

This paper has three new elements. First, we give in the Appendix all analytical formulae 
needed for the calculation. Second, we describe in detail our method for improving the 
full one-loop calculation. As pointed out in [12], in the case of the decay of A — > b\ 62 or 
A — > f\ f<i the decay widths can receive large corrections, especially for large tan (3. This 
makes the perturbation expansion unreliable. In some cases the width can even become 
negative in the on-shell renormalization scheme. We present here a detailed description 
how this problem can be solved by an appropriate definition of the tree-level coupling 
in terms of running fermion masses and running trilinear coupling Af (f = b,r). For 
consistency between the on-shell and the DR running parameters an iteration procedure 
is necessary. Moreover, in addition to the numerical results shown in [12] we present 
here new results for A —> t\t2, A — > bib 2 , A — > f 1 f 2 and the corresponding crossed 
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channels. 



The paper is organized as follows. In section 2 we summarize the tree-level formulae. 
In section 3 the full electroweak corrections are presented for which the explicit analytic 
formulae are given in the Appendices B, C, D. In section 4 a detailed description of 
the method to improve the one-loop calculation is presented. In section 5 the numerical 
treatment as well as numerical results are shown. Section 6 summarizes our conclusions. 



2 Tree— level result 

The sfermion mixing is described by the sfermion mass matrix in the left-right basis 
Cfi> /a), and in the mass basis (fa, / 2 ), / = t, b or f , 

m] L a f?f \ t / m |, 



* \ atmt m~ / ^ ' \ m~ 

\ J J JR / \ .12 

where R ia is a 2 x 2 rotation matrix with rotation angle 9j, which relates the mass 

eigenstates fa, i = 1,2, (m^ < m^) to the gauge eigenstates f a , a = L,R, by fa = R{ a f a 
and 

mj L = M 2 {Q L} + (I 3 f L -e f sin 2 9 w ) cos2/5mJ + m) , (2) 
mj R = Mffj q gy + ef sm 2 9 w cos 2(3 m z + m 2 , (3) 

a f = Af - fi (t&n f3)~ 2I f L . (4) 

Ma, Mi, Mq, Mg and Mg are soft SUSY breaking masses, Af is the trilinear scalar 
coupling parameter, \i the higgsino mass parameter, tan f3 = ^ is the ratio of the vacuum 
expectation values of the two neutral Higgs doublet states [2, 3], Ij L denotes the third 
component of the weak isospin of the fermion f,ej the electric charge in terms of the 
elementary charge eo, and 9 W is the Weinberg angle. 

The mass eigenvalues and the mixing angle in terms of primary parameters are 
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m h, 2 = 2 + m fR T v (m L ~ m V 2 + 4a / m / ) ' (5) 

cosOj = = ~ a/ 2 7/ = = (O<0,<7r), (6) 

and the trilinear breaking parameter Af can be written as 

rrifAf = m 2 LR + mf fi (tan (3)~ 2I i (7) 

with m\ R = (m| - m\) sm9 f ~cos9 f ~. 

At tree-level the decay width of A — > /i/ 2 is given by 

Nl K(m 2 A0 , m 2 ; , m 2 ; ) 

10 7T T72^o 



with K(x,y,z) = \J (x — y — z) 2 — 4yz and the colour factor Nq = 3 for squarks and 

= 1 for sleptons, respectively. G{j 3 denotes the A°—f*-fj coupling as given in Ap- 
pendix A. 



3 Full Electroweak Corrections 

The full one-loop corrected decay width is given by 



(9) 



16 7rm^ 



|G{ 2 3| 2 + 2Re(G{ 2 3-AG{ 23 



The (UV finite) corrections AG{ 23 consist of the vertex corrections SG{^ (Fig. 12), wave- 

f(c) 

function corrections and the coupling counter term corrections <5G{ 23 owing to the shift 
from the bare to the on-shell values, 

AG{ 23 = 5G(^ + 5G{^ + 6G{® . (10) 



The renormalization procedure with the fixings of the counter terms is given in [12]. The 

explicit formulae of the vertex corrections 8G{ 23 as well as the various contributions to 

/(«) 

123 ) 



the wave-function corrections 8G^ X 



8GS = 2^ 
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5Z[ 1 + 5Z( 2 + 5Z^\g{ 23 , (11) 



can be found in the Appendices B and C. For the explicit formulae of the self-energies 

needed for the calculation of the counter term correction, 5G{^ (see eq. (23) of [12]) , we 
refer to the Appendix D. 

Due to the diagrams with photon exchange we also have to consider real photon emission 
corrections to cancel the infrared divergences (Fig. 12). Therefore the corrected (UV- 
and IR-convergent) decay width is 

pcorr^O _^ jjj = T{A ^ jjj + p^O ^ ^ (12) 

Throughout the paper we use the SUSY invariant dimensional reduction (DR) as regu- 
larization scheme. For convenience we perform the calculation in the 't Hooft-Feynman 
gauge, £ = 1. 



4 Improvement of One— loop Corrections 

It has already pointed out in [12] that in the case of bottom squarks and tau sleptons, 
especially for large tan /3 the corrections to the decay widths A — > b\ 6 2 and A — > f\ f 2 
can be very large in the on-shell renormalization scheme. If the corrections are negative, 
the one-loop corrected width can even become negative and therefore unphysical. Hence 
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the perturbation expansion around the on-shell tree-level is no longer reliable. It has 
been shown in [13, 14] that, in the case of the decays into bottom squarks, the source 
of these large corrections are mainly the counter terms for mj, and the trilinear coupling 
A b , in particular the SUSY-QCD corrections. However, despite the absence of strong 
interactions for the decay into tau sleptons, the corrections become extremely large. We 
show that this problem can be solved by defining an appropriate tree-level in terms of 
running values for rrif and Af. The expansion around this new tree-level then no longer 
suffers from bad convergence. 



Correction to mt, 

First we review the improvement of the perturbation expansion by using DR running 
bottom quark masses, following [14, 15, 16]. 

If the Yukawa coupling h b is given at tree-level in terms of the pole mass nib, the 
one-loop corrections become very large due to gluon and gluino exchange contributions 
to the counter term dm^. The large counter term caused by the gluon loop is absorbed 
by using SM 2-loop renormalization group equations in the MS scheme [14, 15, 16]. Thus 
we obtain the SM running bottom rh b {Q)su'- 



rh b (Q) 



MS / ™b{Q) SM 



SM 



m b( m &)sM 



™&(™&)SM 



(13) 



The ratio (mb(Q)^/nib(mb)^j can be expressed as 



rh b {Q) 



MS 
SM 



c 5 (ar(m b )/ir) 

c 6 (a^(Q)M c 5 (aW(m t )/n) 



c 6 (ai 2) (m t )/7r) c 5 (ai 2) (m b ) / ir) 



(m b < Q <m t ) 

(Q > m) , 



where we have used the functions 

1 9 

23 
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— x) (1 + 1.175a;) 
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r 



Cq{X) ■ 

and the 2-loop RGEs for a s [16], 

12vr 



a 



(2) 



(Q) 



(33-2n / )ln^- 



6(153 



(m b < Q <m t ) 
(Q > m t ) , 



(33 



2n f y 



In 



(14) 



A- 



with n f = 5 or 6 for m& < Q < m t or Q > m t , respectively. For the SM DR running 
bottom mass at the scale Q = nib we use the MS equation 



™fe( m &)sM 



1 + 3 



+ A' 



7T 



7T 



(15) 
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with K q = 12 A and then convert to DR using one-loop running a s (Q): 

~ (n \ I m fe(<5)sM \ - / nMS a s(Q) 

m fe (g)sM = — — -= m 6 (m 6 )^ — m fc (16) 

In the MSSM, for large tan /3 the counter term to mj can be very large due to the gluino- 
mediated graph [13, 17, 18]. Here we absorb the gluino contribution as well as the sizeable 
contributions from neutralino and chargino loops and the remaining electroweak self- 
energies into the Higgs-sfermion-sfermion tree-level coupling. In such a way we obtain 
the full DR running bottom quark mass 

^(Q)mssm = rhb(Q)sM + 5m b (Q) . (17) 

The explicit form of the electroweak contribution to the counter term 8m b (Q) is given in 
Appendix D.4. 

Correction to A^^ 

The second source of a very large correction (in the on-shell scheme) are the counter 
terms for the trilinear coupling A bjT , (see eq. 7), 

SAbT = ^ - H^k^l + Sntm/3 + fiStaa/3. (18) 
m b>T m byT m b:T 

Again, the big bottom mass correction Sm b contributes to SA b , but also the counter term 
of the left-right mixing elements of the sfermion mass matrix, Sm 2 LR , gives a very large 
correction for higher values of tan/3. In particular, in the case of the decay into staus, 
this is the main source for the bad convergence of the tree-level expansion. As in the case 
of the large correction to m b we redefine the Higgs-sfermion-sfermion tree-level coupling 
in terms of DR running A& ^m^o) . Because of the fact that the counter terms 5A b<T 
(for large tan (3) can become several orders of magnitude larger than the on-shell A b)T we 
use A biT (rriAo) as input [14]. In order to be consistent we have to perform an iteration 
procedure to get all the correct running and on-shell masses, mixing angles and other 
parameters. This procedure is described below. 

5 Method of improvement 

In this section we will explain in detail how we can improve the perturbation calculation 
for the sbottom and stau case by using DR running values for m b and A bjT in the Higgs- 
sfermion-sfermion tree-level couplings. Since we take DR running values for A b and A T as 
input and all other parameters on-shell we will have to pay attention to the sbottom and 
stau sector in order to get consistently all needed running and on-shell masses, mixing 
angles and other parameters. Here we adopt the procedure developed in [14] and also 
extend it to the electroweak case. 
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5.1 Calculation of running and on— shell parameters 



Stop sector: 

We start our calculation in the stop sector. Because all input parameters in the stop 
sector are on-shell we obtain the on-shell masses m f ~ i; m f ~ 2 and the stop mixing angle Q~ t 
by diagonalizing the stop mass matrix in the t^-tn basis, see chapter 2. The running stop 
masses m^. and mixing angle # f ~ are calculated at the scale Q = Qi= \l m i x m t 2 by adding 
the appropriate counter terms to the on-shell values in 

m|(Q f ) = m| + 5ml, (19) 
5ml = »l4(m|), (20) 

h{Qt) = 0t + S9i. (21) 

The electroweak parts of the sfermion self-energies 11^ (mi) are given in Appendix D.6 and 

the SUSY-QCD contributions, U^ SY ~ QCD (m 2 f ) are given in eqs. (25)-(27) in [8]. Here 

and in the following all running parameters X(Q) are related to their on-shell values X 
by X(Q) = X + 5X, with 5X being the full one-loop counter term — also including the 
SUSY-QCD parts. The mixing angle is fixed by [19] 

60 f = \ (SZ{ 2 - 8ZQ = Re (nf 2 (m|) + uUm\) 1 • (22] 

h h 



For DR running m t we use the formulae from section 4 with the obvious substitutions 
m& — > m t and K q = 10.9 for the top-case. Next we evaluate the running parameters 
Mq(Q) and Mq(Q) by inserting the running values ml(Q),9 i (Q) J m t (Q) M ssM, rh z (Q) = 
m z + 5m z , $(Q) = (3 + 5(3 and 6 W = 6 W - - ^) into the equations 

Mj = m| cos 2 9 f + m 2 i2 sin 2 9- t - m\ - m 2 cos 2(3 (lf L - e t sin 2 6> w ) , (23) 
Mq = m| sin 2 Q~ t + m| cos 2 6~ t - m\ - m\ cos 2(3 e t sin 2 6 W . (24) 
For the running value of A t we use (see eq. (7)) 

A t = (m| - m|) + A cot $ , (25) 

where we have taken running fi(Q) = \i + (5A) 22 [20, 21]. 



Sbottom sector: 

In the sbottom sector we have given all parameters on-shell except the parameter for 
the trilinear coupling, A b (Q), which is running. First we calculate m fe (Q^)MSSM from 
eq. (17) at the scale Q~ b = ^/m^m^ . From the stop sector we already know the run- 
ning values of Mq,ta.n(3 and [i. Then we diagonalize the sbottom mass matrix using 

m b(Q{)MSSM, Mq, tan (3, fi and on-shell Mj>, which is near its running value Mj>, to ob- 
tain the starting values for and The on-shell sbottom masses m~ b . and the mixing 
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angle 6 b are calculated from their running values by subtracting the appropriate counter 
terms, i. e. m~ = rh? (Q) — 5m? , 9 b = 9~ b (Q) — S9 b . Now we can compute the running 
value for M^. Using the relation (see eq. (3)) 

Mq = mj- sin 2 9 b + cos 2 9 b — m\ — m 2 cos 2(5 sin 2 6 W (26) 

we get Mq = (M? + 8M? f' 2 « M 6 + 5M?/(2M^) with 

= (5^771^ sin 2 6^ + ^tti 2 ^ cos 2 6 b + — m|J sin 29 b S9 b — 2m^ 5m;, 

—bm\ cos 2/3 e?, sin 2 9 W + 2m 2 sin 2/3 5(3 e& sin 2 6* w — m 2 cos 2/3 e& 5sin 2 9 W 

(27) 

and = ^(Q^mssm — f^b- Because the parameters involved in these calculations are 
very entangled, e. g. depends on the Sm b , which themselves depend on M^, we have 
to perform an iteration procedure. 

Iteration procedure: 

Here we will describe in detail the procedure how we obtain all necessary on-shell and 
running parameters. For convenience, we shortly denote all masses, parameters, couplings 
etc. for a certain n > 1 in the iteration by X^ n \ As starting values X^ we take on- 
shell masses and parameters (except which is running) and the couplings derived from 
these quantities. The only exceptions are the standard model running fermion masses 

shortly stands for the full DR running fermion masses, ?tij(<5)mssm- 
The single steps of the iteration procedure are the following: 

1. The running stop masses and the stop mixing angle are calculated as explained 
above by m? (n) = m| + 5m 2 ^ n \x^) and Of = 6- t + SOf (X^). 

2. m[ n) =m ttSM + Smi n) (X^) 

3. mi n) = m z + 8mt\x( n -V) and 

sin 2 6$ = sin 2 9 W + 5sm 2 6$ with 5sin 2 6$ = - cos 2 6 W (^- —) (X^) 

\ m w m z J 

4. The running value of tan/3, tan/3^ = tan/3 + 5tan/3^ ?1 - ) , 

with 5tan/3^ = Iran^o^™" 1 )) tan/3 [221. 

m z sm2/3 

5. /i (n) = n + 5/i (n) with SfiW = 5X 2 2(X {n -V). 

6. The soft SUSY breaking masses are calculated from m^\ 9~ , m t (Qt)^ n \ 
m ( z \ sin 2 ^ and tan/3^. 

7. We compute the running A t by using running values in eq. (25): 

A (n) = f^Un) 2(„)X Sin29^ + ^ ^ ^ 

2 Wit 
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8. In the sbottom sector we obtain 8m b from the running values already calculated 
in steps 1.-7., like m~ n \ or rhf, and the remaining masses, couplings etc. from 

Af(n-1)_ 

9. m[ n) = m 6j sM + ^[ n) • 

10. We receive the running sbottom masses, m~ , and the mixing angle, 9^\ by solving 

Oj 

the mass eigenvalue problem with the running values of , , m£ n \ A b , (x^ 

and tan (3^ n \ 



11. The on-shell sbottom masses m 2 ^ = rru — Srrit} ' \Q%^) at the scale - 
mt ] mt\ and 9^ = # n) ~ ^ ■ 

bi 02 b b b 

12. 5Ml (n) = 5m 2 ^ n) sin 2 B~ h + 5m 2 ^ n) cos 2 9 b + (m^-m^J sin 29 b 59^ - 2m b (m [ b n) - m b 

-Sm 2 z (n) cos 2(3 e b sin 2 9 w +2m 2 z sin 2/3 5(3^ e b sin 2 w -m| cos 2(3 e b Ssin 2 6%\ 
(Remember that the values without a hat (~) are on-shell ones!) 

13. M^M fi + -^. 

D 2 M 6 

14. In the sneutrino sector we calculate the running sneutrino mass 
2(n) - -l T + 5m 2 >\x^) and M 2{n) = m 2 > ] - \ cos2/3», see also eq (23). 



m,~. = m 



15. In the stau sector the values for running ml n ),m~ etc. are calculated like in the 
steps 8-13 in the sbottom sector with the evident substitution b — > f for the corre- 
sponding parameters and Mq — > Mi, — > Mg. 

16. All couplings are recalculated with the new running parameters — > X n . 
The iteration starts with n — 1 and ends, when certain parameters are calculated precisely 



enough for a given accuracy, i. e. 



^ _ a-( n ) 



< £ for £ = {m b , Mfy,m T , Mg}. For e we 



choose £ = 10 8 . We have checked the consistency of this procedure by computing the 
on-shell Mjj and running Mq from the sbottom sector by using 

Mp = mg sin 2 6^ + m~ 2 cos 2 9 b — m 2 — m 2 cos 2(3 e b sin 2 9 W , (28) 
M? = rh~ hi cos 2 ^ + mf 2 sin 2 9 b -m 2 b - m 2 cos 2/3 (/ b 3i - e fe sin 2 W ) , (29) 

which are equal (up to higher order corrections) to the on-shell input Mfy and running 
Mq from the stop sector. 

For easier reading the single steps of the iteration procedure of the stop and sbottom 
sector are depicted in the flowchart in Fig. 1. 
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Figure 1: Simplified flowchart for the iteration procedure. For details see section 5.1. 
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5.2 Numerical results 



In the following numerical examples, we take for the standard model parameters m z = 

91.1876 GeV, m w = 80.423 GeV, sin 2 W = 1 - m 2 w /m 2 z , a{m z ) = 1/127.934, m t = 

174.3 GeV, m& = 4.7 GeV, m T = 1.8 GeV and ^m u ,md,rn e ,rn c ,rn s ,'m iJ ) = {4,8,0.511, 

1300, 200, 106} MeV for 1 st and 2 nd generation fermions. M' is fixed by the gaugino 

5 

unification relation M' = - tan 2 9 W M, therefore the gluino mass is related to M by 

m § — {a s {rrig) / a) sm 2 6 W M . In order to reduce the number of parameters in the input 
parameter set, we assume Mq = Mq 3 = fM^ = yfM^ = M Ls = = Mq i2 = 
Mr-, = Mn = Mj = Mfp for the first, second and third generation soft SUSY 

(-'1,2 in, 2 -^1,2 Cj\,2 ' ° 

breaking masses as well as A = A t = = A T for all (s)fermion generations, if not stated 
otherwise. 

stop— case: 

In Fig. 2 we show the tree-level and the corrected widths to A — > tit 2 for tan (3 = 15 
and {m A o, A, M, Mq} = {700,-500,120,300} GeV as a function of the higgsino mass 
parameter [i. The electroweak corrections are almost constant about —7%. 
At /j, ~ —242 GeV one can identify the pseudo-threshold coming from t 2 t x% 




-400 -350 -300 -250 -200 

fi [GeV] 

Figure 2: Tree-level (dotted line), full electroweak corrected decay width (dashed line) 
and full one-loop (electroweak and SUSY-QCD) corrected width (solid line) of A — > £1^2 
as a function of /i. 

Fig. 3 shows the tree-level, the full electroweak and the full one-loop corrected (elec- 
troweak and SUSY-QCD) decay width of A — > t\t2 as a function of the lighter stop 
mass, m^, where Mq is varied from 200 to 450 GeV. As input parameters we choose 
{m A o,/i,A,M} = {900,250,300,120} GeV and tan/3 = 7. Again, in a large region of 
the parameter space the electroweak corrections are comparable to the SUSY-QCD ones. 
The pseudo-threshold at pa 304 GeV originates from £2 - > ^ X3 i n the wave-function 
correction. 

In Fig. 4 the dependence of the crossed channel decay width, T(t 2 — > tiA°), as a 
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Figure 3: Tree-level (dotted line), full electroweak corrected decay width (dashed line) 
and full one-loop (electroweak and SUSY-QCD) corrected width (solid line) of A — > t±t2 
function of mi . 



function of tan/3 is given. We see that the electroweak corrections have different sign 
compared to the SUSY-QCD ones and go up to 10%. As input parameters we have 
chosen {m A o, /i, A, M, Mq} = {170,500,-390,250,350} GeV as well as M 0s = 450 GeV 
to get an acceptable stop mass splitting. 




5 10 15 20 25 30 

tan (3 



Figure 4: tan /^-dependence of the tree-level (dotted line), full electroweak corrected 
(dashed line) and full one-loop corrected (solid line) decay widths of t 2 — > t\A°. 

Fig. 5 shows the decay width T(t 2 — > tiA°) as a function of mf x , varying Mq 3 from 200 
to 460 GeV. To get a larger mass splitting for the top squarks, we relax the conditions 
for 3rd generation squarks and take {M^, M^} = {500,300} GeV. All other SUSY 
breaking masses are fixed at 300 GeV. For the remaining input parameters we choose 
{m A o,/i,A} = {120,-400,-350} GeV and tan/3 = 7. 



12 




250 300 350 400 
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Figure 5: mf -dependence of the tree-level (dotted line), full electroweak corrected 
(dashed line) and full one-loop corrected (solid line) decay widths of t% — > t\A°. 



sbottom— case: 

In Fig. 6 we show two kinds of perturbation expansion for r(v4° — > bib 2 ) with {m^o, /i, 
A, M, Mq] = {800, -300, -500, 200, 300} GeV: First we show the tree-level width, given 
in terms of on-shell input parameters (dotted line). The dashed and dash-dot-dotted 
line correspond to the on-shell electroweak and full (electroweak plus SUSY-QCD) one- 
loop width, respectively. For both corrections one can clearly see the invalidity of the 
on-shell perturbation expansion, which leads to an improper negative decay width. The 
second way of perturbation expansion is given by the dash-dotted and the solid line 
which correspond to the improved tree-level and improved full one-loop decay width, 




30 35 



Figure 6: Two kinds of perturbation expansion: the dotted line corresponds to the on- 
shell tree-level width, the dashed and dash-dot-dotted line correspond to electroweak 
SUSY-QCD on-shell one-loop width, respectively. The dash-dotted line corresponds to 
improved the tree-level and the solid line to the (full) improved one-loop width. 
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respectively. Here we take the same input parameters as in the first case but with running 
A = -500 GeV. 

In Fig. 7 we show the decay width T(A° — > b\ b 2 ) as a function of DR running A for 
{m^o, /i, M, Mq} = {800,-300,300,300} GeV and tan/? = 30. The dotted line corre- 
sponds to the improved tree-level, the dashed line corresponds to the improved SUSY- 
QCD one-loop width and the solid line shows the full improved one-loop width. For 
negative A the electroweak corrections decrease the decay width by ~ 20%, whereas for 
positive A the SUSY-QCD corrections almost vanish and the electroweak ones go up to 
30%. 
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Figure 7: Perturbation expansion around the improved tree-level decay width (dotted 
line) of T(A° — > b\ b 2 ) as a function of the trilinear coupling A. Dashed and solid lines 
correspond to the improved SUSY-QCD and full improved one-loop width, respectively. 
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Figure 8: tan /3-dependence of Y(p2 — > b\A°) for two kinds of perturbation expansion. 
The dotted and dash-dot-dotted lines corresponds to on-shell tree-level and full one- 
loop width, respectively, the dash-dotted line corresponds to improved tree-level and the 
solid line shows the full improved one-loop width. 
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Fig. 8 shows the behaviour of the decay width r(&2 —> biA°) for large tan/3. As in 
Fig. 6 two kinds of perturbation expansion are given. The dotted and dash-dot-dotted 
lines correspond to the tree-level and full one-loop decay widths in the pure on-shell 
scheme. For large tan (3 one can clearly see the invalidity of the perturbation series, 
leading to a negative decay width. In the second case we show the expansion around 
the tree-level decay width, given in terms of running Aj, and m&. The dash-dotted 
line corresponds to the improved tree-level and the solid one to the one-loop width. 
Up to tan/3 ~ 30 the corrections stay relatively small which indicates that already the 
(improved) tree-level is a good approximation for T(b 2 — > b\A°). As input parameters 
we take the values {m A o, //, A, M, Mq} = {150,-220,500,200,300} GeV and M &3 = 
500 GeV for kinematical reasons. 

In Fig. 9 the A decay into two staus is given as a function of tan (3. Despite 
the absence of SUSY-QCD corrections the perturbation expansion around the on-shell 
tree-level (dotted line) leads to an improper negative decay width (dashed line) com- 
ing from large 0(h%) corrections. As input parameters we take {m^o, /i, A, M, Mq} = 
{800,400, —500, 120,300} GeV. The dash-dotted line corresponds to the improved tree- 
level and the solid line shows the improved one-loop width for the same input parameters 
as above and running A = —500 GeV. 
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Figure 9: On-shell tree-level (dotted line) and full electroweak on-shell corrected decay 
width (dashed line) of A — > f\f% as a function of tan/3. The dash-dotted and solid lines 
correspond to improved tree-level and full improved one-loop decay widths. 

Fig. 10 shows the decay width of T(A° — > fi^) as a function of the mass of the 
decaying Higgs boson A for the improved perturbation expansion. The dotted and 
the solid lines correspond to the (improved) tree-level and full one-loop widths, respec- 
tively. In the whole region of the parameter space shown the (electroweak) corrections 
decrease the on-shell width by 15%. As input parameters we choose {//, A, M, Mq} = 
{-450, -500, 120, 260} GeV and tan/? = 7. 
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Figure 10: m A o -dependence of the improved tree-level (dotted line) and full one-loop 
corrected (solid line) decay widths of T(A° — > fif^). 

In Fig. 11 we show the A dependence of r(f 2 — > fiA°) in the improved case. For 
negative values of A the corrections increase the on-shell width by ~ 20% whereas for 
positive values of A the corrections are negative and go up to 15%. The input parameters 
are taken as follows: {m A o,jj, M, Mq} = {150,400,300,300} GeV and = 500 GeV 
for an acceptable stau mass splitting. For tan/3 we take the value 30. 




Figure 11: A-dependence of the improved tree-level (dotted line) and improved one-loop 
decay width (solid line) for f 2 — > fiA°. 
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6 Conclusions 



We have calculated the full electroweak one-loop corrections to the decay widths A — > 
fi f'2 and J2 — > f\A° in the on-shell scheme. We have presented all formulae required 
for the computation. It has been necessary to renormalize almost all parameters of the 
MSSM. We have also included the SUSY-QCD corrections which were calculated in [8]. 
For the decay into bottom squarks and tau sleptons for large tan (3 an improvement of the 
on-shell perturbation expansion is necessary. We have worked out an iterative method to 
improve the one-loop calculation. Thereby, the tree-level coupling is redefined in terms 
of DR running masses and running Af. We find that the corrections are significant and 
in a wide range of the parameter space comparable to the SUSY-QCD corrections. 
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A Notation and Couplings 

For the neutral and charged Higgs fields we use the notation H k = {h ,H°,A°,G }, 
Hi = {H + ,G + ,H~,G-} and H k = (Hi)* = {H~ , G~ , H+ , G + }. t/t stands for an 
up-type (s)fermion and b/b for a down-type one. Following [2, 3] the Higgs-Sfermion- 

Sfermion couplings for neutral Higgs bosons, G{j k , can be written as 

G{ jk = G(H»Jff 3 ) = [RfGi R>k {Rf) T ],.. (30) 

The 3rd generation left-right couplings G^ LRk for up- and down-type sfermions are 

G l f -V2h t m t c a + g z m z (If L -e t s 2 w )s a+l 3 —^(A t c a + fis a ) 

LR ' 1 I -^(A t c a + fiSa) -V2h t m t c a + g z m z e t s 2 w s a+ p 



G 



( V2h b m b s a + g z m z (Il L -e h s 2 w )s a+ p %( A b s a + fj.c a ) 

h h t A „ , \ . fK-i „, „ i „ ™ „ „2 



LR ' 1 I ^|(A b s a + yUc Q ) \/2h b m b s a + g z m z e b s 2 w s a+f} 

Glr2 = g lri with a — > a — 7r/2 , 
G* Li?3 = —V2h t 

G LR,3 = ~V2h b 



-fAtC^ + zis^ 



2 



-±(A b sp + ncp 

| (A 6 S/3 + /i Cp) 



GL,4 = GL.3 with/3^/3-7T/2, 

where we have used the abbreviations s x = sinx, c x = cosx and s w = sin# w . a denotes 
the mixing angle of the {h°, if °}- system, and h t and h b are the Yukawa couplings 

= -7= -■ (31) 



w sin P \p2m w cos /3 

The couplings of charged Higgs bosons to two sfermions are given by (I = 1, 2) 

Gg' = G(Hff*f^j =GQ = (RfGg^Rff)^, (32) 
/ /if,m;, sin P + /itTOf cos P — fl^Bp gin 2/3 WA&sin/3 + /xcos/3) \ 

G LR1 = ; ( 33 ) 

\ hf(Af cos P -\- fi sin /3) h t m b cos P + h b m t sin P J 

f h b m b sin P + h t m t cos P — sin 2 P h t (A t cos/5 + /x sin/3) \ , - ,t 
I /ib(^4fe sin/3 + /J cos/3) h t m b cos P + h b m t sin P J ^ lr,ij 

(34) 

G{1 2 = with - - | . (35) 
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/' denotes the isospin partner of the fermion /, i. e. t' = b, b\ = ti etc. Note that only 
the angle (3 explicitly given in the matrices above has to be substituted; the dependence 
of (3 in the Yukawa couplings has to remain the same. 



The H®Hif*fj interaction is given by 



kl I e ij 



H k H l Ji J j , 



(36) 



with 



c kl 



kl 



sin 2 a 
— \ sin 2a 



/ cos 2 a 
| sin 2a 





- \ sin 2a 
cos 2 a 







\ 



sin 2 /? 

■\ sin 2(3 cos 2 /? J 



-|sin 2(3 



\ sin 2a 
sin 2 a 










\ 






cos 2 /? 
\ sin 2(3 sin 2 /? J 



1 sin 2(3 



(37) 



(38) 



2c 2 , 



{if - e f S w) R il R jl + e f S w R i2R S j2 



?/ uf 



f nf 



(39) 



For the H^H t f*fj interaction 
1 _ 

£ = 



J2 [h) d f kl (R{ 2 R f j2 + r{[rQ+ 9 2 (di - 4)4] //,; //, /;/ 
/ 



J ■ 



we use the coupling matrices 

/ sin 2 (3 -\ sin 2(3 
,1 — \ sin 2(3 cos 2 (3 



ki 



d 



ki 













V 













sin 2 (3 - 1 sin 2(3 
-±sin2/? cos 2 /? J 



( cos 2 (3 ±sin2/? 
\ sin 2(3 sin 2 /? 








\ 



V 



o 
o 



o 








cos 2 /? 
| sin 2(3 sin 2 /? J 



1 sin 2/? 



(40) 



(41) 



(42) 



J i j 



f 



2cl 



If cos 29 w - e f s 2 w ) RiiRjx + e f s z w R{ 2 R J j2 



f of 



?/ T)f 



For the Higgs-fermion-fermion couplings the interaction Lagrangian reads 

2 4 2 

£ = E 4 A*// + E4 fifc/T 5 / + E [^ + * {ViPr+ViPl) b + h.c. 



fc=i 



fc=3 



1=1 



(43) 



(44) 
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with the couplings 



- 9 



rat cos a 
2mw sin 3 



2mw sin B 



-^cosa, 



^sina, 



rrif, sin a 
" 2m\Y cos 8 

mfc cos a 



^sina, 



5 



2mn; cos 8 



1 9 



y\ 



2mw 

2m w 
mt cot 8 



i^sin/3, 



V* = 9—2 



mt 



rn\\- 



h t sin j3, 



112 



a 2mw 

2m w 
m b tan 8 



^sin/3, 
- z-^|cos/3, 



hb cos /3 . 



(45) 



The interaction Lagrangian for Higgs bosons and gauginos is given by 



^ - -o Z^ H kXl^lmkXm- }^ H kXl flmkJuX 
Z k=l Z fc=3 



rO-Oj-,0 ~0 -5 



m 



-</ £ if fe ° ^ + (f+,P r + P £ ) x+ + ty £ Pf fc ° %f (f&P* + F+ ik P L ) x) 



k=l 
2 



fc=3 



-g £ [#+ ^ + (P* P* + F&P L ) xl + h.c. 



(46) 



k=l 



with 



Imk 



F + 

ijk 



2 

dk 

+T 

1 



ZuZ m 2 + Z m iZi2 — tan0 w {Z\^Z m \ + Z m ^Z\\ 



F° mlk , (47) 



(48) 



and 



F R 
1 ilk 



ilk 



VaZ lA + —^( z i2 + Z n tan6> w )^ 2 



■efc+2 



i3 



71' 



Z t2 + Z a tan0w)^ 



(49) 



U, V and Z are rotation matrices which diagonalize the chargino and neutralino mass 
matrices, dk and take the values 

d k = { — cos a, — sin a, cos fl, sin/3} , = {— sin a, cos a, — sin (3, cos /?} . 
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The coupling of the vector boson Z° to two sfermions, £ = —i g z z(j Z® f* <9 M fj with 
g z = g/ cos9 w , is given by the matrix 

z{j = C[ R^R^ + C{ t R{ 2 Rj2 
with C[ = Ij L — tjS 2 w and C R = —ejS 2 w . 

For the interaction of a vector boson with two gauginos we use the couplings 
Ofj = Zi 2 Vji — -^Zi4Yj2 , Of- = ZiiUj\ + -^ZisUj2 , 

O'if = -VaVji-^V^ + Sijsl, 
{ f = — UnUji — \Ui2Uj2 + SijS^ , 

= — 7;Z i3 Zj3 + -Z^Zj^ = — . 

The interaction Lagrangian of the chargino-sfermion-fermion couplings is given by 
C = i (l*.P R + klP L )xfb i + b (f ij P R + l&P L )xT t ii 

+ xj ($Pl + klP R )tb* + xf (fi 3 P L + ^P R )bt* (50) 
with the coupling matrices 

1% = -gV jx R\ x + h t V j2 R\ 2l l\ = -gU n Ri + h b U j2 R\ 2 , 

(51) 

k% = hU^Rl, k% = h t V 32 R\ x . 

For the neutralino-sfermion-fermion couplings the Lagrangian reads 

C = f (a{ k P R + b{ k P L ^j x°Ji + X k (o&l + H k P R ) ff* (52) 
with the coupling matrices 

a ik = hjZ kx R f i2 + gf[ k R f n , b ik = h s Z kx R s iX + gf Rk R{ 2 (53) 

and 

fL = V2((e s -lf)t^6 w Z kl + lfZ k2 ) , f f Rk = -V2e f t^6 w Z kl . (54) 
x takes the values {3, 4} for {down, up}-type case, respectively. 

B Vertex corrections 

Here we give the explicit form of the electroweak contributions to the vertex corrections 
which are depicted in Fig. 12. For SUSY-QCD contributions we refer to [8]. 

+ 5G{^ + + SG&rt + 8G(ir Z) + (55) 
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The single contributions correspond to the diagrams with three scalar particles (5G{^' H ^ 
and 5G[z^ HH ^ , three fermions f^Crf^j'* and dG^l'^^, one vector particle (SGl^^ 
or two scalar particles (SGi^' and 5G{^' H J in the loop. SG^^™^ denotes the 
correction due to the mixing of A and Z° and 5G[^ AG) is the Higgs mixing transition 
A°-G°. 

As shown in [12] we can sum up the A°Z° and A°G° transition amplitudes which leads to 

5G{£ AZ) + 5G{il> AG) = -—U AZ (<o)G{ 24 . (56) 

The explicit form of the A°-Z° self-energy, Il^^m^o), is given in app. D.l. The vertex 
corrections from the exchange of one Higgs and two sfermions are 



5G{^' Hff) = -JT-T5 E E G L3GLfe^.C , o(m|,m^o,m|,m^o,mJ m ,m 

V^l m,n=l k=l 



~ 77^ E E G Cn2,G f imk G "nk C (™>}. , m\ , m| , 171%+ , 771^ , m| ) 
\ / m,n=l k=l r<-n\ 

(57) 

with the standard two-point function Co [23] for which we follow the conventions of [24]. 
The graph with 2 Higgs particles and one sfermion in the loop leads to 

^ G 123 = -J^y^T^ E (EE G imk G mjl A k,l-2 + E E G imk G mjl A,k-2 ) x 



m=l \ k=l 1=3 k=3 1=1 



s~i ( 2 2 2 2 2 2 

Co m.-, m A o,rrif ., m? ,m H o,m H o 

\ J i J J J TIT, £ 



— — /J L m w E ( G im2 C'o (m| , m^o ,mj.,m~ L ,m 2 H +, m 2 G+ ) 



m=l 



/-<//' /"<//' n ( 2 2 2 2 2 2 

~ (j im2 U jml °0 ^ m j. > m A j rn f.j' rn fr i m G+ ' m H+ 



with 



.4 



ki 



I — cos 2(3 sin(a + (3) — sin 2/3 sin(a + (3) \ 
\ cos 2/9 cos(a + /?) sin 2/3 cos(a + (3) ) 

For the gaugino exchange contributions we get 

- 1 4 - 

8G(il' xff) = — — - E F(m 2 f.,m 2 A0 , mj m f , m f ; s f 3 , -s(, b{ k , a{ k , a f jk , b f jk 
y^j k=i 

+ J^y X! F \ m % m A^ m % m xii m f'' m f'i s 3 si', k{ k , l{ k , l f jk , k f jk ) , 



k=l 
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1 4 - - 

k,i=i 

(58) 

where F(. . .) shortly stands for 
F(ml m 2 , m 2 , M , M u M 2 ; ^, g*, g{, g*, g%) = (h 1 M 1 + h 2 M 2 )B (m 2 , M 2 , M|) 
+ (M4o+/iiMi)S (m?, M 2 , M 2 ) + (/loA/c + ^M,)^™ 2 , M 2 , M 2 2 ) 
+ [2(^f^f+^^)M MiM 2 + /i Mo(M 2 + M 2 -m 2 ) + ^(M^+Ml-m 2 ) 
+ /i 2 M 2 (M 2 + M 2 -m 2 ) ] C (m 2 , m 2 , m 2 , M 2 , M 2 , M 2 2 ) (59) 

with the abbreviations h = (g% gf gf + godi^) , h i = (9o9i92 + S'W #2 ) and h 2 = 
(9o9i92 + #0 <?f <?2 ) • For up-type sfermions F kl3 = F^ l3 and for down-type sfermions 
chargino indices are interchanged, F^ 3 = Fj^ 3 . 

We split the irreducible vertex graphs with one vector particle in the loop into the single 
contributions of the photon, the Z-boson and the W^-boson, 

8G^ V) = 5G{^ + 8G^ Z) + 6G{% W) . (60) 

In order to regularize the infrared divergences we introduce a photon mass A. Thus we 
have 

SG i23 7) = r^y ( e oe/) 2 G{ 23 v(m 2 f.,m 2 A0 , m|, A 2 , m|, m|) , 
1 2 

5G£' Z) = — — g\ Gins 4m 4j V (m| , m\ , m| , m 2 , m 2 ^ , mjj 

V m,n=l 
2 2 

~ 71^ G Uk Z L Rxk{a-P)V (m\ , m 2 ~ , m| , m| , m 2 H o , m 2 ^ 

v 47 v Z k,m=l 
i 2 2 

+ Ta^V ~f Z ° fimk z LjRik(a-P) V (mj , m\ , m 2 ^ , m 2 , m 2 ^ , m 2 H o 

V 471 "/ Z fc.m=l 



fc,m=l 
.2 2 



Arx f>,W) _ 1 9 V- r /' d/d/d/' d/' T// rr ,2 _2 
0Lt 123 _ 77 \ 2 IT" ^mn3 - n 'il- ft jl- n 'ml- n 'nl \ fi ' A 

l 47F J Z m,n=l 



,2 _2 _2 _2 _2 



J 3 Jm Jr. 



„2 2 . _ 

9 /^f//' d/ d/' T/Y 2 2 2 2 2 2 

(47rp 2^2 ^ G omi R n R mi V [m A o,m f -^m f ^m w ,m H+ ,m^ 



* g 2 2 



(4 
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where we have used the vector vertex function 



V (ml m 2 Q , mi Ml Ml M% 



-Bo (ml Ml Ml) +B (ml Ml M 2 ) + B (ml Ml M 2 ) 



+ (-2ml+m\+ml-Ml + Ml + Ml)C (mimlmlMlMlMl) (62) 



and the rotation matrix R 



ki- 



R k i(<f>) 

z{a can be found in Appendix A. 



cos <p sin i 
— sin d> cos < 



ki 



For the vertex graphs with 2 sfermions in the loop we obtain 



SG 



/(«■//) 

123 



(Aw) 2 
1 

"(4^ 



m,n=l 
m,n=l 



f 

nm'S 



nm3 * 



p/ I p/ 

^ijmn ' -^mnij ly CX 1 Hnmj 



{2lf-e f )e f si)R% mn + eW w R, 



f 

mjin 



B (l 



m A o,m? ,m 



fn 



2 J2 p/fl 



tjmn 



(JVg + 1) 



+ e f) e f s w Nc(yR{j mn +R mni ^+R{ nm: j+R : 
x £? ( % , m? , ml 



/ 

mjin 



1 



C ~ ~ ~ £ 

Nckfhf ^ G^ m3 (Rijnm + Rjimn) ^0 (j 



m A o , m ; , m = 



( 47r ) m ,n=l 

For various products of sfermion rotation matrices we have introduced the short forms 



(63) 



R 



h 
ijkl 



p/ p/ p/ p/ 
rt il- rt jl- rt A;l- n ';i ) 

?/ p/ p/ p/ 



p/i? p/ p/ p/ p/ 

-'V/fci ~~ rL i2 rl j2 rL k2 rL l2 ' 



i? 



/ p/ p/ p/ p/ 

ijkl — rL il rL jl rL k2 rL l2 5 

//f p/ p/ p/ p/ 



(64) 



p/ pj p/ p/ 

Hjkl — rL il ri j2 rL kl rL l2 ■ 

Note that the last term in eq. (63) originates from the mixing of 2 squarks and 2 sleptons, 
where / denotes the 'family partner' of the fermion / with the same isospin and from the 
same generation, i.e. t = v T or fj = &j. The diagrams with one Higgs boson and one 
sfermion in the loop lead to 

4 2 



5G 



f(v,Hf) 



123 



(An) 2 



k=3 m=l 



G Lk ( h) c{ k 5 mj + g 2 (c b 3k - c* 3fc ) e f mj j B (m 2 A0 , m 2 H o , m\ 



(4tt) 



k,rn=l 



(h 2 -g 2 /2)cos 2 f3-(h 2 -g 2 /2)sm 2 {3 
(h 2 f + Zip — g 2 ) sin /3 cos /? 



p/' p/ 



+hfhf> 5 k2 R^R^ 



BQ[m A0 ,m H+ ,mr, 



with /if = {/i<,0} and /ij = {/i^, /i r } for the decay into {squarks, sleptons}, respectively. 
The Higgs-Sfermion coupling matrices c kl and e{j can be found in Appendix A. 
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Figure 12: Vertex and photon emission diagrams relevant to the calculation of the virtual 
electro weak corrections to the decay width A — > f\f 2 - 
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C Diagonal Wave— function corrections 



For the diagonal wave-function renormalization constants we use the conventional on- 
shell renormalization conditions which lead to 



5Z 33 



SZL 



f 



(65) 



where the dot in Uu(k 2 ) denotes the derivative with respect to k 2 . In the following we 
list the single contributions of the wave-function corrections (Fig. 13 and 14). 

C.l Higgs part 



(Att) 2 



N c (4) B (m 2 A0 ,m 2 ,m 2 f ) + m 2 A0 B (m 2 AO ,m 2 f ,m 2 ) 



(4tt) 2 
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(72) 
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.2 Sfermion part 
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Figure 13: Diagonal Higgs self-energies 




Figure 14: Diagonal sfermion self-energies 
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D Self— energies and counter terms 



Here we give the explicit form of the self-energies needed for the computation of various 
counter terms for the one-loop width A — > fif%. 

D.l AZ— mixing 

The scalar-vector mixing self-energy, U A z(k 2 ), is defined by the two-point function 




M = -ik^U AZ {k 2 )tl{k) 



The single contributions from the particles x are denoted by the superscript x in H AZ . 
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Figure 15: AZ-mixing self-energies 



D.2 W + self— energies 

For the calculation of the mass counter term of a gauge boson V (V — W ± , Z°), 5m v = 
3?Ilyy(my), we need the transverse part of the vector self-energy Uy V (k 2 ) from 
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Here /-[■ and fi denote up- and down-type (s)fermions of all three generations, respectively. 
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Figure 16: W + self-energies 
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D.3 Z Q self— energies 

Accordingly to eq. (83) the mass counter term contributions to the gauge boson Z° are: 
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D.4 Fermion self— energies 

In our notation, the fermion self-energy is denned by 




M = iu(k)U(k)u(k) 



with 



n(A;) = f^p L u L (k)+ jiP R u R (k) + u SL (k)p L + u SR (k)p R . 

Thus the counter term for quarks and leptons is given by 
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Note that for quarks and leptons (contrary to charginos), the left- and right-handed 
scalar parts of H(k) are equal, U SL (k) = U SR (k). The single contributions to 5mj are as 
follows: 
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Figure 18: Fermion self-energies 



D.5 Chargino self— energies 

Since the higgsino mass parameter /i is fixed in the chargino sector, the counter term S/i 
reads [20, 21] 
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with the chargino self-energies IT^ = IL^m 2 ^). U and V are two real 2x2 matrices which 
diagonalize the chargino mass matrix, 

UXV T = M D = 



m~+ 

Xi 

m~+ 

A2 



The single left- and right-handed parts of IT^ can be found by comparing the coefficients 
accordingly to eq. (101). 
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Higgs/gaugino contribution: 
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Figure 19: Chargino self-energies 



D.6 Sfermion self— energies 

For the fixing of the sfermion mixing angle 8j we need the off-diagonal elements of the 
sfermion self-energies, = IT^m^ ). In the following, Y[j R denotes the weak hyper- 
charge, y[/ R = 2(lj L ^ R — ef). The short forms for various products of sfermion rotation 
matrices can be found in B. Additionally, we use the abbreviation R^§ = RliRj 1 R k2 Ri2- 
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The contributions from first and second generation sfermions, F m , are given by 



W = WU' - F ') > i'//" = n& r (/' - /") 
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(125) 



where the sub-/superscript F denotes values belonging to first and second generation 
scalar fermions with same isospin as / (e. g. F% = {ui,C\} for the stop case, ...), F' 
sfermions with different isospin etc. 
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Figure 20: Off-diagonal Sfermion self-energies 
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